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Abstract. We derive the eikonal equation of light wavefront in the presence of Lorentz invariance violation 
(LIV) from the photon sector of the standard model extension (SME). The results obtained from the 
equations of E and B fields respectively are the same. This guarantees the self-consistency of our derivation. 
We adopt a simple case with only one non-zero LIV parameter as an illustration, horn which we find two 
points. One is that, in analogy with Hamilton- Jacobi equation, from the eikonal equation, we can derive 
dispersion relations which are compatible with results obtained from other approaches. The other is that, 
the wavefront velocity is the same as the group velocity, as well as the energy flow velocity. If further we 
define the signal velocity v s as the front velocity, there always exists a mode with v s > 1, hence causality 
is violated classically. Thus our method might be useful in the analysis of Lorentz violation in QED in 
terms of classical causality . 
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1 Introduction 

The constant speed of light in vacuum is one of the basic assumptions of special relativity, and it serves as a crucial 
ingredient of Lorentz symmetry, which is a cornerstone of modern physics. However, currently there is a revival of 
interest in the possibility of varying speed of light, or in other words, a tiny deviation from exact Lorentz invariance, 
in the search of so-called quantum gravitational phenomena. 

As the typical quantum gravitational scale, i.e., the Planck scale -Epianck = a/cS/Gc 2 , is practically unattainable 
from conventional accelerator experiments, people turn to search for phenomenologically accessible effects from quan- 
tum gravity [I] at relatively low energy scales. One of them is the possibility of Lorentz invariance violation (LIV). 
The possibility that quantum gravity may leave a tiny imprint of LIV at relatively low energies was observed by many 
authors from various approaches to quantum gravity. These include string field theory, where the tachyon field may 
induce an instability to the naive Lorentz invariant vacuum and translate it into the potential of a tensor field. As a 
consequence, the tensor field acquires a vacuum expectation value [2] and breaks Lorentz invariance. Later, by incorpo- 
rating various LIV coefficients into the standard model, the theory was developed into an effective field theory (EFT), 
called standard model extension (SME) [3 . Other approaches include spin network calculation in loop gravity [3], de- 
formed special relativity [5] , foamy structure of space-time [5] , noncommutative field theory [7] , emergent gravity [5] , 
and the recently suggested Horava-Lifschitz gravity [9] . All of them suggest that tiny LIV may be a signature of new 
physics. 

Fortunately, LIV not only resides in theoretical considerations, it also becomes experimentally testable, and has 
already been tested to very high accuracy in various sectors of the standard model (for a review, see Ref. [TU]). From 
an experimental viewpoint, people have already been able to severely constrain the linear order LIV correction to the 
photon group velocity, which is suppressed by the ratio of experimental energy E to a large mass scale, e.g., -E/Mpi anclt 
(see e.g. Refs. [TUII1II1II1IIS]). 

From the EFT viewpoint, conventionally non-renormalizable operators can be naturally suppressed by a large mass 
scale, where new physics might come in. So the mentioned tiny linear correction may come from dimension-5 opera- 
tors |16j . However, if one assumes CPT symmetry and the anisotropic scaling between space and time [5] (also termed 
weighted power counting [T7]), or alternatively imposes CPT invariance and the supersymmetric constraints [TS], one 
can naturally expect that the leading order correction comes from dimension-6 operators, thus can evade current ex- 
perimental constraints [TH] ■ On the other side, there is no clear theoretical reason to pin down the values of dimension-3 
and -4 LIV operators, though they have already been extensively studied and found to suffer severe constraints from 
experiments [10P20j. Hence it is still valuable to consider them both theoretically and experimentally, and it inspires 
us to revisit the renormalizable LIV operators [3]. In this paper, we mainly focus on the CPT-even part of the photon 
sector in the framework of SME. 

The paper is organized as follows. In section 2, we give a brief review on the renormalizable photon sector of SME. 
With some ad hoc assumptions for simplification, we derive the modified Maxwell equations as our starting point. 
In section 3, following the derivation provided by Fock [3T], we obtain the eikonal equation of the modified Maxwell 
equations implicitly. In section 4, by adopting a much simpler case with only one non-zero LIV parameter as an 
illustration, we give the eikonal equation explicitly. From the derivation, we find it is self-consistent, because one can 
obtain the same result from equations relevant to either an E or B field, respectively. It is also consistent with other 
approaches, as the dispersion relations obtained from our procedure and those from others are the same. Meanwhile, 
in this simple case, we find that the group velocity equals the wavefront velocity, as well as the energy flow velocity. 
As the derivation of eikonal equation is obviously not restricted to the simple model, the method may be useful in the 
classical causality analysis in LIV extension of quantum electrodynamics. Section 5 summarizes the results obtained 
in this paper and gives a further discussion on LIV topics. 



2 A brief review of the photon part of SME 

SME, including in principle all possible LIV terms in the framework of EFT [3], is a well-motivated testable approach 
to LIV physics, and it has been extensively tested in various sectors of the standard model (TUj. Within SME, the 
renormalizable photon sector is one of the most tested parts [20,22 j. Its Lagrangian reads 



photon = —^F^F^ - ~{k F ) KXliV F KX F^ + -(k AF ) K e KX ^A x F^, (1) 

where the coefficients k F and k AF characterize the violation of Lorentz symmetry. From Eq. ([1}, the equation of 
motion is deduced as 

d a F lia + (k F )^ aM d a F^ + (k AF ) a e^ 7 F^ = 0. (2) 
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With the 3+1 decomposition, we define 

(k DE y k = -2(k F )°i° k , (k HB f = l(k F y kmn e^ k e lmn , 

(k DB y k = -{k HE ) k * = ^(k F )°^ mn e kmn , (3) 

where the Latin indices run over the three spatial coordinates, from 1 to 3. After subtracting the trace part of the first 
two matrices, we define further the following traceless matrices which are frequently used in this paper, 

I3e ij = (k DE y 3 - a6V, fa* = {k HB ) ij + aft, j ij = (k DB ) ij = -(k HE ) l ° , (4) 

where 

a = -tr(k DE ) = --tr(fcHs). (5) 

Note that the {k^^p-y tensor has the same symmetry as the Riemann tensor, R^apt [31123]. Due to the double 
tracelessness of k F and the Bianchi identities (k E ) ^[vpa] = 0, we can get Eq. (O and trfy] = 0, respectively. The first 
two matrices in Eq. (jj) can be easily shown as symmetric, i.e., {fi E ) T — (3 E and (/3b) T = Pb- Since the (Auf)^ term 
may cause theoretical instabilities [3] and has suffered severe astrophysical constraints [53], we abandon it from now 
on. Moreover, as 7 y mixes E and B fields and may cause further complications in calculation, we simply set 7 y = 0. 
Further, we also set a — as an ad hoc assumption. Under these simplifications, there are only 10 parameters left, 
residing in the traceless symmetric matrices f3 E and j3 E - Finally, we can write down the Lagrangian in terms of E and 
B as 

£ P hoton = \{E 2 - B 2 ) + X - (((3 E y k EJE k - (f3 B y k BiB k ) , (6) 
where the electric field E and the magnetic field B are defined conventionally as 

B i = -±j* k F jk , E l = —d t A i — S7i(f> = F i0 . (7) 

From Eq. ([5]), we can deduce the modified Maxwell equations, 

V-E + di(/3 E ) ij E j = 0, -E i -(/3 E ) i JEJ+6 ijk ajB k + e ijk aj(/3 B ) kl B 1 = 0, (8) 

and the equations 

V-B = 0, VxE + B = (9) 

come from the 3+1 decomposition of Bianchi identity <9 M *F A " y = 0, where *F txu = ^e^^ 1 Fp-y- The above four equations 
are our starting points of derivation in the next section. 

Note that Eq. © has the analogy of electrodynamics in a homogeneous anisotropic media [3] if we regard 
as the inverse of magnetic permeability and (f3 E ) l: > as the dielectric constant, i.e., as defined in Ref. [33], D 1 = 
(<5 1J + /3E 1J )Ej, H 1 = (£ 1J + /3B 1J )Bj. From this analogy, we expect that the differential equation satisfied by the light 
wavefront in LIV vacuo should have a similar form as the eikonal equation in an anisotropic medium. Actually, we 
will see that this is indeed the case in section 4. 

Now we turn to the discussion of the energy-momentum tensor, which is useful for the discussion of causality 
in section 4. The conventional procedure to obtain a symmetric energy-momentum tensor is through the so-called 
Bclinfante tensor [5S] 



where 



However, this does not work in the presence of LIV, as pointed out by Colladay and Kostelecky [3J. 
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In fact, from Eq. (jlOp . we can get the energy-momentum tensor corresponding to the Lagrangian in Eq. (JT|) without 
the CPT odd term, 

&»> = -[F» a F v a + {k F ) a ^ s F v s F al3 ] - if v C 

= -[F pa F» a + {k F ) a ^F u s F aP ] + \v^[(kF) a lS F a0 F^ + F a ^F aP ]. (12) 

In the above derivation, we have used the matrix representation of the homogeneous Lorentz algebra for a covariant 
vector field, 

(D;=v/-V). (13) 

We find that due to the presence of the second asymmetric term, — {k F ) a ^ p& 'F v s F a p, one can no longer obtain a 
symmetric energy-momentum tensor. This may be one feature of LIV theories. One consequence is that the definition 
of the conserved 4- momentum density is Op , instead of pQ . Actually, with direct calculation it is easily checked that 
d p pp — 0, hence d(J d 3 x0 o,1 )/dt — (for details, see Appendix A). In contrast, this is not valid for 0^°. The spatial 
part of the latter is denned as the generalized Poynting vector [3]- For the convenience of our discussions in section 4, 
we give those components of the energy-momentum tensor explicitly in terms of E, B fields, 

6)00 = \{E 2 + B 2 ) ~ {k F f j0k E' 3 E k + -(k F ) ijkl e ijm e kln B m B n 

= i(£ 2 + B 2 + E ] {(i E y k E k + B rn ((3 B ) mn B n ), (14) 

m = -F^F^ - 2(k F ) 0l0k F l jF ok - {k F f ikl F l 3 F k i 

= e ljk E j B k + e ijk (/3 E ) jl E l B k , (15) 
0i o = _pij F o. _ (k F y ikl F° 3 F kl - 2(k F y j0l F° j F 0l 

= r* 1 E'ir + e l]k E ] {(i B ) kl B l , (16) 

where we have used the reparametrization of Eqs. Q and and the simplifications corresponding to the Lagrangian 
in Eq. ([6]). From Eqs. (IT5l) and fll"6|) . it is apparent that due to the loss of Lorentz invariance, characterized by Pe and 
j3 B , generally the energy-momentum tensor is not symmetric, i.e., 0^ v — V>1 = \(k F ) al3llS F^ — {k F ) a ^ vS F^ s ]F a p 7^ 0. 

It has been well known that, in the Lorentz invariant theory, the energy-momentum tensor should be symmetric 
due to the following reasons [26] : 

(1) When it couples with gravity in curved space-time, the symmetrical Einstein tensor G^ v and metric tensor 
automatically imply a symmetric 0^ v due to the Einstein equation; 

(2) The current conservation of M ppv = \{x v Q p ^ - x p pv ) of the angular momentum tensor = J d 3 xM 0pu also 
requires a symmetric pv . 

However, in our case, LIV undermines these two reasons; 

(1) It is necessary to take LIV effects into account in gravity side for a consistent theory, however in this case, it is 
the local Lorentz invariance involved. Moreover, when coupled with gravity, LIV parameters should be promoted 
to dynamical fields instead of being constants [37] , and the compatible framework to incorporate LIV in general 
relativity is the Riemann-Cartan geometry |28j . where the extended Einstein tensor is no longer symmetric. The 
antisymmetric part of the Einstein tensor and the energy-momentum tensor is consistently linked through field 
equations [28] . Thus no inconsistency arises when an asymmetric p,v couples with gravity in a spontaneous LIV 
theory [28]. 

(2) Due to the breaking of Lorentz invariance, the generator of the Lorentz group, i.e., the angular momentum tensor, 
is no longer conserved in general, thus no constraints are imposed here. If one assumes that only boost invariance 
is broken, i.e., rotational invariance is still preserved. It means that only a and k° AF are non-zero in Eq. ([IJ. Then 
in this case, one can prove <9 y = 0^ 1 , which is required by the conservation of the spatial part of the angular 
momentum tensor, i.e., dE 1 ^ jdt = 0. However, in the presence of LIV, it is only valid in a particular reference 
frames, e.g., CMB reference, as the transformation from one inertial frame, in which the LIV photon equation is 
isotropic, cannot remain isotropic in other inertial frames [29j . Here we present the proof in Appendix A. 



3 A derivation of eikonal equation 

In order to derive the eikonal equation of the modified electrodynamics, it is necessary to review how this is achieved 
in Lorentz invariant Maxwell theory. We follow the derivation of Fock [3T] . The basic principle is that the fields on the 
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wavefront, which is a space-like 2-dimension surface, must be noncontinuous in the normal direction of the surface. 
The reason is that the surface is an interface of two regions, one already influenced by the field perturbation and the 
other remaining intact. Thus the time derivatives of the fields on the wavefront should be singular. 

Taking a right-moving plane wave as an example, the wavefront is an infinite plane perpendicular to its constant 
wavevector k at any instant. The left region of the plane is a region with non-vanishing electromagnetic fields, while 
the right side is a region with vanishing fields. This conveys the fact that light propagates with a finite speed, thus 
at any instant it cannot propagate further to the right of the wavefront plane. Since an electromagnetic field evolves 
dynamically, it is more helpful to analyze from a 4-dimensional point of view, rather than from the view at some fixed 
instant in 3-spatial dimensions. 

For any instant t, the wavefront in the 3-spatial dimension can be written as 

t = f(x,y,z). (17) 

This can also be viewed in a dynamical way as an evolving 2-dimensional space-like surface f(x, y, z) — t = 0. When 
we write Eq. p7[) in an implicit form, 

w(t; x, y, z) = 0, (18) 

the wavefront becomes a slice of 3-dimensional hyper-surface in flat space-time in the 4-dimensional point of view. 
So the 2-dimensional space- like wavefront is just a time slice of the null-surface (light-cone), and the non-continuity 
of fields on the wavefront is nothing more than the (classical) causality requirement. To be more specific, the fields 
inside the light-cone must be continuous since they are causally connected. Using Maxwell equations, one can obtain 
from the fields on any given surface in 3-spatial dimensions inside the light-cone, the fields on an infinitesimally close 
surface. However, this cannot be done for fields lying on the surface outside the light-cone. So the instant time slice of 
the hyper-surface w(t; x, y, z) — becomes an interface between causally connected and non-connected regions, hence 
the derivatives of fields on the interface become singular. We use the above principle to derive the so-called eikonal 
equation with the same logic as that in the Lorentz invariant theory. 
The Lorentz invariant eikonal equation reads 

(V/0r,y,z)) 2 =n 2 , (19) 
where in vacuum, n equals 1 in our unit (i.e. c = 1). In a more explicit Lorentz invariant form, it reads 

(^) 2 -(V-) 2 =0. (20) 



We show in Appendix B how to get Eq. (|20p from Eq. (Q3 

Now we turn to derive eikonal equation in the presence of LIV. Any field u(x,y,z,t) on a 3-dimensional hyper- 
surface in the 4-dimensional Minkowski space is 

u(t, x, y, z) = u(f(x, y, z); x, y, z) = u (x, y, z). (21) 

Therefore, we have 

du du du df 

_ = . _| i_ (22) 

dx 1 dx l dt dx l ' 

where x % with Latin index i running over 1 to 3, represents x, y, z respectively. Thus the spatial derivatives of the 
electric and magnetic fields on the hyper-surface are 

From the linearity of Eq. (f2"B"| , we have 

V x E = V x E + Vf x E. (24) 

Then by using the second equation of the Bianchi identity we get 

V x E = -B + Vf x E. (25) 

Similarly, one can obtain from Maxwell Eqs. (JSJ) and (J9j) the set of equations below 

V ■ ((1 + /3e) • Eo) = Vf ■ ((1 + (3e) • E), VxE = -B + Vf xE; (26) 
V-B = Vf-B, Vx ((1 + /3 B )-B ) = ((l+/3 E )-E) + Vf x ((l+feJ-B). (27) 
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Then by multiplying the second equations of ([25]) and (|2"T)) respectively with V/, we have 

V/ • V x E = -Vf • B, (28) 
V/ • V x ((1 + p B ) ■ Bo) = Vf • ((1 + E ) • E). (29) 

Comparing the equations above with the first equations of ([26)) and ((27|) , we find the following relations 

V/ ■ V x E + V • B = 0, (30) 
V/ • V X ((1 + B ) ■ Bo - V • ((1 + E ) • E ) = 0. (31) 

When f(x, y, z) is a constant, i.e., by specifying a particular instant, Eqs. (1301) and (|31l) become the first equations 
of the original field equations ((8]) and ([9]) respectively. This means that we need to specify proper initial conditions. 

In order to extract more information, we multiply the curl of Eq by the tensor (1 + 0b), then calculate its cross 
product with V/. By utilizing the second equations of (12^)) and (|27|) successively, we get 

V/ x [(1 + B ) • V x E ] = Vf x [(1 + B ) • (Vf x E - B)] 

= (1 + 0b) • E — V x [(1 + B ) • B ] + Vf x [(1 + B ) • (Vf x E)]. (32) 

After arranging it in such a way that all time derivatives of E are on one side, and Eq and Bq on the other side, we 
have 

(1 + 0b) • E + Vf x [(1 + B ) • (Vf x E)] = Vf x [(1 + B ) • V x E ] + V x [(1 + B ) ■ B ]. (33) 

Now it is apparent that the right-hand side of (|33[) includes the wavefront function / and the electromagnetic fields 
Eq and Bq with values on the wavefront, while the left-hand side includes only E and /. With given wavefront function 
/ and the fields Eq and Bq, E on the interface must be singular, as we have already argued. Otherwise one can obtain 
the fields on an infinitesimal close surface outside the light-cone from field equations, which obviously violates the 
causality requirements. Thus, if we view the left-hand side of Eq. (|33p as a matrix equation, 

M e ■ E = (1 + E ) • E + Vf X [(1 + B ) • (Vf x E)], (34) 

where M e is a tensor/matrix, then the determinant of M e must be zero, i.e., Det(M e ) = 0. Otherwise, in principle, 
one would get a non-singular E by solving Eq. ((33]) . 

Similarly, we can obtain the corresponding equation of the B field, 

B + Vf x [(1 + 0b)- 1 • {Vf x ((1 + B ) • B)}] = Vf x {(1 + 0b)- 1 • (Vf x [(1 + B ) • Bo])} - V x E , (35) 

where we have presumed the existence of the inverse of (1 + 0b), due to the fact that all of the 10 parameters, (0b) 1 - 7 
and (PbY 3 , must be very tiny compared to 1 to meet the stringent experimental constraints. 
One can also define 

M b ■ B = B + Vf x {(1 + 0E)- 1 ■ (vf x [(1 + B ) • B]) } . (36) 

And for the same reason, Det(M{,) = 0. Actually, we find that generally Det(Mb) cx Det(M e ), hence the two equations, 
Det(Mfc) = and Det(M e ) = 0, give the same differential equations of f(x,y,z). The former proportionality (in some 
cases, it even becomes an equality) holds as expected, as the differential equation of / should be unique. This will be 
shown in detail in Appendix C. In the presence of LIV, Det(Af e ) = actually gives the eikonal equation implicitly. 
We will show it explicitly in a simple case in the following section. 



4 A case study and discussions 

In this section, we give explicitly the eikonal equation for a simple case, together with more detailed discussions. For 
this purpose, Eqs. (fM)) and (|5o| are our starting points. 

First we try to extract the tensors/matrices M e and Mb from Eqs. (|34p and (|36[) respectively. By rewriting Eq. (IMl) 
in component form, 

{M e ym = (1 + /3 E fW + e i J k f i (l + B ) kl e lmn f m E n , (37) 
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where fa = df/dx 1 , one can read 

(M„)« - [1 - (V/) 2 ]^' + fifj + (0 E ) ij ~ e mk ^ ml f n f m ((3 B ) kl . (38) 
Similarly, from Eq. f|36[) . one can obtain 

(M 6 )« = ^ - ^"'^"/n/^" - e lsfc e" m 7 s / m ^ M (/3 s )" J , (39) 

where W™- 7 = [(1 + Pe)^ 1 ] 13 ■ Our task now is to calculate the determinants of M e and Mb, then one can obtain the 
eikonal equation satisfied by f(x,y,z). 

As E, B fields are simply the 3+1 decomposition of the electromagnetic field strength , and our working 
hypothesis relies only on the analysis of causality, we can reasonably expect that the eikonal equations obtained from 
the equations of E and B should be the same. 

The explicit expressions of M e and Mb in the matrix form are rather tedious, not to mention the calculation of 
their determinants. As an illustration, we need not to consider the full expression with all 10 parameters non-zero. We 
place the discussions of the full expression in Appendix C, where we show that the eikonal equations obtained from 
Eq. ([55)1 and Eq. (j3"9")l are indeed the same. Instead, we assume here that only (Pb) 12 = (Pb) 21 = cr ^ to simplify 
our discussions on M e and Mb- By virtue of this simplification, we can subsequently get and solve the eikonal equation 
explicitly. 

However, here we derive it by restarting derivations from the simplified field equations, 

V-E = 0, e ijk 9 j B k + e ijk 5 j (/3 B ) kl B 1 -E i = 0, (40) 
VB = 0, e ijk 6> j E k + B i = 0, (41) 

where (Pb) H = <r(S k 5 2 + S^5\). Now the equations corresponding to Eqs. ([55)1 and (|33p arc 

B - (Vf) 2 [(l + fa) • B] + Vf [Vf • [fa • B)] = Vf x {V x [(1 + fa) • B ]} - V x E - Vf(V • B ), (42) 
[1 - (V/) 2 ]E + Vf x {fa • (Vf x E)] = Vf x [(1 + fa) ■ (Vf x E )] + V x [(1 + fa) • B ] - (V • E )Vf. (43) 

From (|4"5]l we obtain 

(Me) - [1 - (V/) 2 ] I 1 I + a I ./':? -hh I • (44) 




By direct calculation, we have 

Det(M e ) = [1 - (V/) 2 ] {[1 - (V/) 2 ] 2 + 2/r/ 2 a[l - (V/) 2 ] - fl<J 2 (Vf) 2 } . (45) 
Similarly for the B field, from (|4"2")l we have 

/ioo\ / hh A 2 -(v/) 2 o\ 

(M b ) = [1 - (V/) 2 ] 1 + a /| - (V/) 2 hh , (46) 

\ooi/ \ hh hh oj 

and 

Det(M 6 ) = [1 - (V/) 2 ]{[1 - (V/) 2 + ahhf - (fi + fl)(fl + flV)- (47) 

One can easily verify that Det(Af e ) = Det(-Mft). Generally, the equation Det(M e ) = has three solutions. One is the 
conventional (V/) 2 = 1, and if written in terms of w(x, y, z; t), it is Eq. ([2"U]) . We will see below that this corresponds 
to the conventional dispersion relation p 2 — 0, which accompanies modified Lorentz-violating dispersion relations in 
many models, see e.g., Ref. [12] , The other two solutions, 

1 - (V/) 2 + ahh ± ^(/i 2 + /!)(/! + /!) = 0, (48) 

manifest LIV explicitly, which is characterized by the non-zero LIV parameter a. Treating space and time on the same 
footing, i.e., rewriting Eq. (|48p in terms of function w(t,x,y, z), we can get 



A 2 - (Vw) 2 + a(w lW2 ± y/(wl + wl){wl + wl)) = 0, (49) 
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where Wi = dw/dx 1 . When a = 0, Eq. (14T)1) reduces to the Lorentz invariant case, i.e., Eq. (|2T))) . 

In order to solve the eikonal equation, we use the Lorentz invariant case (|20p as an illustration. Choosing the 
positive sign for convenience, we obtain a first order partial differential equation, 

dw 



<)t +V(Vw) 2 = 0. (50) 
Then from the well-known analogy between Hamilton- Jacobi equation and geometric optics |21j . one can identify 



w(x,y, z;t) as the Hamilton action S(q,P;t), (Vw) 2 = y X)i=i ^1 as the Hamiltonian H(q, fj|;t), and i«i as the 
momentum Pj = dS /dq l . So the question is transformed into solving the Hamilton equation, 



We note that X^iO*'*) 2 = 1> w hich can be written in an appropriate form, dr 2 = dt 2 — (dx 1 ) 2 = 0, characterizing 
the null geodesic of photons. If Eq. (|50[) is solved in the momentum space, we get the well-known dispersion relation 
p 2 = for massless particles. 

Next we use the same analogy to solve Eq. (PRO) . In the presence of LIV, the formal Hamilton- Jacobi equation now 
reads 



dw 
~dt 



+ J (Vw) 2 - a(w!W 2 ±yj(w{+ w 2 ){wl + wl)) = 0, (52) 



with 



H(iVi) = J (Vw) 2 - a{w x w 2 ±\J(w\ + w 2 )(wl + w 2 )). 



(53) 



Since H does not contain the canonical variables which conjugate to momentum Wi, it automatically implies that the 
momentum fc, = Wi is conserved. From the absence of explicit dependence on time, we have dS/dt — dw/dt = —E = 
— fco. Then based on the above observations, we can quickly read the dispersion relation 



k° =k 2 ai^k 2 ± yj (k l2 + P 2 )(fc 22 + fc 32 )), (54) 

which can also be confirmed by solving field equation ^ (for details, see Appendix D). 

Note that the "±" signs in Eq. (jM]) imply that the photon dispersion relation depends on polarization, hence can 
lead to the so-called vacuum birefringence effect [3"ll2~4"]. which has already been used to place severe constraints on 
LIV parameters from astrophysical observations, like CMB [31], radio galaxies [35] and GRB [3TJ1IT5] . 

Back to the Hamilton- Jacobi equation (|52j) . by using the observer rotational invariance, we can rotate the observer 
frame to a particular one with w\ — = p. This simplifies Eq. (|52[) to 



dw 



dt 

For the positive sign in (|55|) . we can solve 



+ sj {Vw) 2 - a(p 2 ± (p 2 + w 2 )) = 0. (55) 



. 1 . 2 ldH pyl — cr .3 dH W3y/1 — a 



x = x = -— = -^= x A = - — = r (56) 

2 dp ^J2p 2 + w 2 dw 3 ^2p 2 + w 2 V ' 

where the factor 1/2 comes from the symmetry w\ = w 2 — p. We can read the group velocity from Eq. (|56[) as 
v g = yE^i^ = \/l — cr, or equivalently dr 2 = (f — a)dt 2 — dx 2 = 0. Accidentally, the phase velocity equals the 

group velocity, i.e., v p — H/\Vw\ = yl — o = v g . 
For the negative sign, we have 

± l ± 2 Igg P ±3 dH w 3 (l + a) 

2 dp ^2p 2 + (I + cr)w 2 ' dw z ^2p 2 + (f + a)^ ' 
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Thus dr 2 = dt 2 



\dx\ 2 + dx2 2 + dx^ 2 /(l + a')] = 0, and the group velocity reads 




2p 2 + (1 + a) 2 w 2 
2p 2 + {l + a)w 2 



1 + 



2p 2 + (1 + a)w 2 



cr(l + <j)w 2 




2|Vw| 




2 ' 



where in the last step we have approximated it to the first order of a. Similarly, the phase velocity is 




1 + 



2p 2 + w 2 




2|Vw| 2 ' 




We see that only in the first order approximation of cr, v p = v gi while considering higher orders, v p ^ v g , contrary to 
the case of positive sign. So in general v p ^ v g in the presence of LIV, and this is easy to understand. As the presence 
of the background tensor Uf breaks the Lorentz invariant vacuum, the vacuum is no longer an isotropic and dispersion 
free medium. In this case, just like the behavior of the electromagnetic wave in the conventional dispersive medium, 
generally v p depends on the photon momentum k, i.e., dv p /dk ^ 0. Consequently, v g — dk° / dk = v p + k dv p /dk ^ v pi 
where k = \k\. 

In this case, it is interesting to define an effective refractive index from the comparison of Eq. (|4"8)) with Eq. (TK)|) . 



by comparing Eq. (|54l) with that of the conventional massive particles. 

Then from Eqs. ((55)) and (f5"9")) . we see that, for a > 0, the "+" sign gives n 2 s > 1 and m 2 ff < 0. The analogy with 
optical media tells us that v p < 1, which is confirmed in our special case where v p = y/1 — a. However, the analogy 
with massive particles by means of introducing an effective mass squared m 2 s breaks down. In the conventional case, 
when we have a negative mass squared, it means that we have expanded the theory at an unappropriate point, a 
point which is not a true vacuum. When we quantize the theory around this unstable vacuum, we get a tachyon 
with v > 1 and hence violate causality. But this analogy breaks down here, as the underlying physical mechanism is 
completely different from that of the false vacuum in a Lorentz invariant theory. For the "— " sign, we have n 2 g < 1 
and m 2 s > 0. From the analogy with optical media, we expect v > 1, while naive analogy with conventional massive 
particles indicates v < 1. We see from the special case above that v p — v g m 1 + aw 2 |Vw|~ 2 /2 > 1, which is again 
consistent with the optical media analogy. For a < 0, one can obtain similar results. Thus we conclude that the analogy 
between LIV vacua with optical medium is more reasonable and helpful. 

In fact, our derivation of the eikonal equation is inspired by and starts from the analogy between the electromagnetic 
wave propagating in LIV vacuo and in anisotropic media. As the definition n 2 s is along the same line of the analogy, 
it is natural to expect that the qualitative results (of velocity) are consistent with direct calculations. However, this is 
not true for the analogy with conventional massive particles. The reason is that here the effective mass squared is no 
longer a free parameter in Lagrangian; instead, it is momentum-dependent. More importantly, now Lorentz invariance 
is broken, thus it is not necessary for the maximum attainable velocity of the effective massive particle to be c = 1. 

Moreover, we see that, just as what happened in the conventional dispersive medium |34j . in our case, the definition 
of light velocity is more involved than that in the conventional Lorentz invariant vacuo. Perhaps it is more complicated, 
for now it is the vacuum itself becoming anisotropic. Though in the dispersive medium, there already exist cases with 
v p > c or v g > c |34j . they do not conflict with the requirement of causality. Since in practice, we can only measure 
directly the velocity of light signals, denoted as v s . The only difficulty in that case is a proper definition of signal 
velocity. Once we have an appropriate definition of v s , it can be shown that v s < 1 [33]. Thus no causal problem 
arises there. Since in that case, our starting point, the wave equation, does indeed obey Lorentz invariance. Only the 
medium in which light propagates, is no longer the Lorentz invariant vacuum, but crystal or nontrivial QED vacuum. 
The medium then singles out a preferred direction, which leads to the so-called soft breaking of Lorentz invariance |35) , 
with a super luminous phase/group velocity. 

However, this is not the situation here. Now it is the basic wave equation, Maxwell equation, which is no longer 
Lorentz invariant. Hence there might be more difficulties in the proper definition of signal velocity and it may cause 
causal problems. For the simple case above, we can define v s as the wavefront velocity Vf. Then since our calculation 
deals directly with electromagnetic wavefront, we find v g — Vf. Alternatively, we can also define energy transport 
velocity v l e — M /0 OQ as the signal velocity. Then by choosing the appropriate ansatz Ai(t,r)± — W±cos[w(t — s ■ 
r/v p ±)] (where "±" sign refers to the two modes in Eq. (|54"]l . respectively) and the Coulomb gauge V • A = 0, we can 




(58) 



Similarly, it is instructive to define an effective mass 




(59) 
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prove \v e \ = v g for both "+" and "— " modes in Eq. (|5"4")l (for details of this issue, see Appendix D). Thus for the simple 
cases corresponding to Eqs. ((56)) and (l57|) . and for either sign of a, we see that there always exists one mode with 
v s > 1, which indeed violates causality. However, note that all the treatments above are purely classical. Since LIV is 
believed to be a quantum gravitational phenomenon, the best way to treat the corresponding causality problem is to 
use quantum field theory. Hence it is more appropriate to discuss the microcausality from the calculation of 2-point 
correlation functions of observables at space-like separations. However, this is beyond the scope of the present paper. 



5 Conclusion 

Inspired by the analogy between the behavior of LIV electrodynamics and electrodynamics in macroscopic anisotropic 
medium, we derived the eikonal equation for Lorentz non-invariant vacuum from the modified Maxwell equations of 
SME. The results obtained from the equations of both E and B fields are the same, and the general case is given in 
Appendix C. This implies that the derivation is self-consistent. Then we use the well-known analogy between geometric 
optics and the Hamilton- Jacobi equation, and find that the solution of the eikonal equation (|48p in the momentum 
space, i.e., Eq. (|54[). turns out to be the modified dispersion relation of photons, which is the same as those obtained 
from other approaches. This fact confirms the consistency of our approach further. 

With the definition of an effective refractive index n c g, we find that the dependence of the velocity on the LIV 
parameter is consistent with the velocity dependence on the refractive index n of ordinary refractive medium. As our 
calculation deals directly with wavefront, the results show that the front velocity equals the group velocity for each 
LIV mode. From the analogy of the LIV vacuum with macroscopic anisotropic medium, we find that the presence 
of LIV makes the proper definition of signal velocity more complicated, especially when considering the existence of 
multi-definition of velocity. Therefore, the treatment of causality becomes much subtle. However, if we define the signal 
velocity as the wavefront velocity, i.e., v s — Vf, naive analysis shows that, for either sign of a, causality is violated 
classically at least for one mode. Of course, to see whether or not tiny LIV might threaten microcausality, one should 
turn to the more trustable and rigorous treatment of quantum field theory to calculate the correlation functions at 
space-like separations. Actually, this was done for the Chern-Simons-like term [37] and it has caused extensive debates 
in the literature (see e.g., Ref. |37|). Complete treatment of microcausality was also achieved for the massive fermion 
sector of SME 3,36, and it was shown that no crucial problems arise for a spontaneous Lorentz invariance breaking 
theory. If Lorentz violation happens intrinsically at high energies rather than spontaneously, as suggested in the so- 
called weighted dimensional model, Anselmi pointed out that it is Bogoliubov's criterion of causality, rather than the 
value of correlation functions, that makes sense [55] . 

As a by-product of this work, we find that the asymmetric energy-momentum tensor 0^ causes no inconsistency 
in principle. By assuming rotational invariance, the conservation of Noether charge, i.e., dSij/dt = 0, requires that 
the spatial components of 0^ v must be symmetric. Through inspection of the corresponding LIV parameters in the 
presence of rotational symmetry, we find that indeed tJ = Jl , though this is valid only in a particular reference 
frames, as already mentioned. While from the definitions of 0^ v and the velocity of energy transport, in our simple case, 
explicit calculation shows that the group velocity also equals the energy transport velocity for each mode. Together 
with v g — Vf, it makes the definition of signal velocity as wavefront velocity more reliable in the classical discussion of 
causality. The last point we would like to stress is that the derivation of the eikonal equation is not restricted to the 
simple model here. In fact, at least for gauge invariant Maxwell equations without E, B mixing terms, it should work 
as well. Thus the method might be helpful in the causality and phenomenological analysis of various LIV extensions 
of Maxwell equations. 
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Appendix A 

This appendix gives the proof that, even for a generic asymmetric energy-momentum tensor, the requirement of 
rotational invariance implies the symmetric relation 0^ — 0^ 1 , which ensures the consistency with the corresponding 
Noether current. Here, the conserved charge corresponding to the Noether current is the angular momentum. 
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For a Lorentz invariant theory, the conservation of generators of Lorenz group is manifested through the Noether 
current, 

M% = i(©^-0>,). (60) 
One can construct the corresponding generators of Lorentz group as 

= J d?xM\ u . (61) 

Then from the current conservation, 

d p M% = 0, (62) 

it is easy to show that dS pv / dt = 0. Hence the six generators of Lorentz group are conserved. Meanwhile, combined 
with the conservation of energy-momentum tensor, i.e., d p O pp = 0, the current conservation implies that, 

o = d p Mw = ^[x^d p pu - x v d p 9 p » + {e pv ^ p - e pp s;)} 

= 1(0^-0^), (63) 

i.e., the energy-momentum tensor is symmetric. However, as mentioned in the main text, generally this is not valid 
in the presence of LIV. 

Here we note that the same reasoning above can also be applied to the rotational symmetry. If we only break boost 
invariance, then the infinitesimal symmetry transformation 

¥(x) -> &{x) ~ l -w^(M l m ^ m (x) (64) 

under Lorentz invariance is replaced with 

!P'(x) f \x) - -w^i^y^ix) (65) 

under rotational invariance. It is nearly the same as ([M]) except that w pv is replaced by (where as before, Latin 
indices i, j run over the three spatial coordinate labels, usually taken as 1, 2, 3; and Greek indices [a, v run over the 
four space-time coordinate labels, 1, 2, 3, 0). 

With the above observation, it is apparent that from the Noether theorem, we can obtain the corresponding Noether 
current M p ^ and Noether charge Eij. With the same reasoning, d p M p% i = also implies <9 U = & 31 . 

The symmetrizability of 0^ is valid for a generic LIV theory with rotational symmetry SO(3) unbroken, like the 
photon Lagrangian with an anisotropic scaling in Ref. [12 . However, we stress again that this symmetric property and 
SO (3) invariance is valid only in a specific inertial reference frames, because LIV implies the existence of a preferred 
direction. While for a specific Lagrangian ([1]) without CPT odd terms, we have 

O pv = - [F pa F» a + {k F ) afipS F v s F aP ] - rf"C, (66) 
guv _ e ^ = _ [^aPMSpv _ (k F ) aPuS F p s ]F aP (67) 



so 



and 



= \F a p[{k F r^ 5 F i s -{k F )^ i5 F^) 

= \Fki[{{k F ) kl]m F\ n - (k F ) kUm F^ m ) + {(k F ) kl ^F l Q - (k F ) kH0 F j )} 

+F k0 [((k F ) kQjl F\ - (k F ) km F\) + ((fc F )^'°n - (k F ) k0i0 F j )]. (68) 
Since rotational invariance requires only a ^ 0, from ([3]) and (U), we have 

(k F ) kljm = -(S kj 8 lm -S km S lj ), (k F ) kho = 0, {k F ) k0i0 = --S ki . (69) 
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Substituting (1BT))) into (16"51) . one can easily find that 

M = ^[(F im Fi m - F jm F\ n ) + (F i0 F{ - F j0 F* )} = 0, (70) 

which is consistent with the requirement of rotational invariance. The same property of %3 can also be checked for 
more general theories with rotational invariance, e.g., the Horava-Lifschitz theory. 

Appendix B 

This appendix shows how to get ([20)) from (fT9|) . 
From (jTTJ) , we have 

dt 

dt = fdx 1 — = (71) 



While from (flgj) . we get 



dt id 

= dw = w i dx i + w t dt => — - = -, (72) 

OX 1 Wt 



where Wi = dwjdx 1 and Wt = dw/dt. Then from (|7ip and (|T2[) . one can easily obtain V,/ = —Wi/wt- So we have 

i = (v/) 2 = ][>,M) 2 = (^) 2 /M 2 - (73) 

i=l 

Appendix C 

In this appendix we show that the eikonal equations we got from ((38|) and (|39|) are the same, and hence are consistent 
with the general arguments in the main text: as the derivation roots in the classical causality analysis, and E, B 
fields are components of the 3 + 1 decomposition of F^ u , hence being connected with each other through the Maxwell 
equations, the eikonal equations derived from the equations of E and B must be the same. 

First, let us focus on the case (/3b) 1 -' = 0. The consistent check of this simpler case will be a little different from the 
case (/3e;) ij ^ . For (Pe) 1 -' — 0, we denote the corresponding matrices of and (|3^|) as M e B and M&s respectively, 
i.e., 

(M efl )« - [1 - (V/) 2 ]^' - e mk e^ ml f n f m (f3 B ) kl , (74) 
(M 6fl )« = [1 - (v/) 2 ]£« - (V/) 2 ^ - fif k {p B ) kj . (75) 

If we assign 




Pb = \h b 5 b 3 , (76) 



the matrices above take the explicit matrix forms as below 

( [1 - (V/) 2 ] + /|(i>4 + i> 5 ) - flH + 2f2f3 b 3 /|l>l - f2fl(>>4 + i>5 ) - /3(*>3/l + t> 2 f 2 ) /l /'3 b 5 + b 2 /? - /2C3/1 + b lH) 

( M eB) = /| bl - /2/1CM + i> 5 ) - /3< i '3/l + *>2/2) [1 - (V/) 2 ] + / 2 (i>4 + b 5 ) - / 2 *>4 + 2/i/ 3 b 2 / 2 /3 b 4 + Ht\ - /l(l>2/2 + *>l/3> 

V + l>2/3 - /2C63/I + &l/3) /2/3 b 4 + 63/? ~ /l< b 2/2 + i>l/ 3 ) [1 - ( V /) 2 ] - (/| b 4 + b 5 ff ) + 2f 2 f ± bj 



Al - (1 + b 4 )(v/) 2 ] + /l(i>4/l + i>l/2 + *>2/3) /l(i>l/l + *>5/2 + i>3/3) - i>i(v/) 2 /l(i>2/l + *>3/2 - (i>4 + b 5 )/ 3 ) - b 2 (v/) 2 \ 

< M 6B> = / 2 (i>4/l + bl/ 2 + i> 2 /3) - b!(v/) 2 [1 - (1 + b 5 )(V/) 2 ] + f 2 C&l fl + &5 f2 + H ft) /2 C&2 /l + i>3 /2 - (*>4 + b 5 )/ 3 ) - b 3 ( V f) 2 

\ ft(&4/l + bl/ 2 + i> 2 / 3 ) - b 2 (V/) 2 /3<i>l/l + 65/2 + f> 3 /3) - b 3( v /) 2 [1 - (1 - i> 4 - i> 5 )('/) 2 ] + ftfe/l + i>3 Si - (b 4 + i> 5 )/ 3 ) ) 

Note that (M e e) T = M e s, i.e., M e s is a symmetric matrix, while M^b is not. So it is not a trivial check that 
Det[Mes] = Det[Mf,s] by direct calculation. Imposing the requirements that their determinants be equal to zero, one 
can obtain the same equations, as expected. Since the results are tedious, we do not present them here. On the other 
hand, if we allow only (/?b) 21 = (Pb) 12 7^ 0, one can easily get back to (|4"5")) and (|4"7} . This can be regarded as another 
consistent check. 

For the case (Pe) 1 ^ 0, as the calculation involves the inverse of (1 + /3_e) u , i.e., W %d = [(1 + Pe)^ 1 ] 1 ^ , the results 
will be more tedious than the previous ones. So we also do not present the details here. We just make two remarks. 
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(1) As the calculation of (|39p involves the matrix W l \ while that of (1551) does not, the determinants of (|38l) and (1551) 
are indeed not the same, contrary to the previous case. However, note that the requirement of the same eikonal 
equations is derived from the vanishing of the determinants of matrices (|38[) and (|39p . This requires only that the 
determinants of these two matrices are proportional to each other, not necessarily being the same. In fact, one can 
show that they differ only by a constant, i.e., 

Det[M e ] = Det[M 6 ] ■ Det[l + P E }. (77) 

So as long as Det[l + /3g] ^ 0, the differential equations obtained from Det[M e ] = and Det[Mf,] = must be the 
same. 

(2) By taking (3e = 0, the solution reduces to those obtained from ([7^)1 or ()75[). except for a missing factor 1 — (V/) 2 . 
This can be traced back to solving field equations, and our calculation is implicitly equivalent to that procedure. 
In solving Maxwell equations, in order to obtain the independent equations of the two physical degrees of freedom, 
we need to choose a gauge, e.g., Coulomb gauge, to eliminate the gauge degrees of freedom (see Appendix D). 
In this process, we eliminate one polarization and leave a constraint, which in turn gives an identity. Eventually, 
we are left with only two independent equations, corresponding to the two physical polarizations. The eliminated 
identity in this process then corresponds to the lacking factor in our method. In our derivation, it is the fifj factor 
in ([55)1 (or the —e mk e^ ml f n f m W factor in (|55)) ') playing the role of gauge constraints to remove 1 — (V/) 2 in (|77|) . 

In conclusion, we see that, as expected, the calculation from either ([55)) or (155)) indeed leads to the same eikonal 
equations, and the derivation is also consistent with solving field equations. Thus it is natural to get from the eikonal 
equations two independent dispersion relations, corresponding to two physical polarizations, see e.g., (|54l) . 



Appendix D 

In this appendix we give an explicit calculation of the energy- momentum flow velocity v\ = l ° /(9 00 . 

Using the ansatz A^(x) = e fJ, (p)ex-p[—i(p°t — p ■ x)) and choosing the Coulomb gauge V ■ A = 0, we get <f> = as a 
special solution from V • E = (i.e., the first equation of 05)) ). Then by substituting Eq. mto Eqs. (HU1) and (|IT)). 
we find that <J4TJ) is satisfied automatically. Then from the left equation of (j40j), we get 

-p 2 Ai + ap 3 (p 2 A 3 - p 3 A 2 ) = 0, 
-p 2 A 2 + ap 3 (piA 3 - p 3 Ax) = 0, 

-(p 2 + 2ap lP2 )A 3 +a(p 1 A 2 + p 2 ^i)p 3 - 0. (78) 
By imposing the Coulomb gauge piAi — 0, the equations above can be reduced to 

V + <rp lP2 a(p\ +pl) \ rAl 



I 2 , 2^ 2 , I { a } = 0- ( 79 ) 

Note that, in this appendix, we do not distinguish the upper and lower indices, e.g., A\ = A 1 . From the existence 
of non-zero solutions of (f79")) , one can easily get the dispersion relation (IM1) . To calculate the energy- momentum flux 
velocity, we use the real component of our previous ansatz instead, i.e., A\(t,r)± = W±cos[w(t — s ■ r/v p ±)}, where s 
denotes the unit vector pointing to the direction of propagation and v p ± denote the two independent phase velocities 
of the two modes in (|M)l . With the help of (f79")) , one can express A in terms of A\, i.e., 

A 2 = -^±^- Al , 4, = _ >^+^i) . (80) 

0-O3+P2) P3 

Substituting these back to ([7]), we can obtain the explicit forms of E and B, where we have assigned 

s = (^cos[0], -=cos[6], Bin[0]). (81) 
v 2 v 2 

Then by substituting the explicit forms of E and B into ([T5)l and (IT3)l . we find that indeed \v l e ± \ — v g ± for each mode. 
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